DTP/99/108 
December 1999 



OS 
On 
OS 



O 

Q 

SO 



X 



The i? ratio in e + e , the determination of a(M|) and 
a possible non-perturbative gluonic contribution | 



A.D. Martin, J. Outhwaite and M.G. Ryskin 
Department of Physics, University of Durham, Durham, DH1 3LE, UK. 



Abstract. We review the determination of the QED coupling at the Z pole, which 
^ is a crucial parameter for electroweak theory. We include recent e + e~ — > hadron data 

£N| , from Novosibirsk and Beijing to re-evaluate a(Mi). We find o(M|) _1 = 128.973 ± 

■ 0.035 or 128.934 ± 0.040 according, respectively, to whether inclusive or exclusive 

£N| ' e + e — — > hadron data are used respectively in the interval 1.4 J~ y/s ^ 2.1 GeV. 

£\J ' The error is mainly due to uncertainties in the data in the low energy region, 

y/s < 2.5 GeV. We find that no advantage is obtained by analytic continuation of the 
• dispersion relation into the complex s plane. We show that the hints of structure for 

' \fs ~ 2.5 GeV may be evidence of a non-perturbative gluonic contribution to R(s). 

23: 

The value of the QED coupling at the Z pole, a(M§), is the poorest known of 
the three parameters (Gf, Mz, a(M^)) which define the standard electroweak model. 
Indeed it is the precision to which we know a(Af§) which limits the accuracy of the 
' indirect prediction of the mass Mh of the (Standard Model) Higgs boson jjj. 

The value of a{Mz) is obtained from 



a- 1 = aiO)- 1 = 137.0359895(61) (1) 
using the relation 

a(s)- 1 = (l - Aa lcp (s) - Aa^ d (s) - Aa^(s)) a~\ (2) 

where the leptonic contribution to the running of the a is known to 3 loops 

Aaiep(Mf) = 314.98 x 10~ 4 . (3) 

From now on we omit the superscript (5) on Aahad and assume that it corresponds 
to five flavours. We will include the contribution of the sixth flavour, Aa top (M|) = 
—0.76 x 10 -4 , at the end. To determine the hadronic contribution we need to evaluate 

as „ f°° R(s')ds' 
Aahad(s) =-—P V (4) 

37T J 4m 2 s'(s' - s) 

I Contribution to the Proceedings of the UK Phenomenology Workshop on Collider Physics, Durham, 
UK, 19-24 September 1999, to be published in J. Phys. G. 
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at s = M|, where R = a(e + e~ — > hadrons)/(j(e + e~ — > ^i + /i~). 

The early determinations (^] of Aahad(Af|) made maximum use of the e + e~ 
measurements of R(s), using the sum of the exclusive hadronic channels (e + e~ — > 
2ir, 3ir, . . . KK, . . .) for yfs < 1.5 Ge V and the inclusive measurement of a(e + e~ — > 
hadrons) at larger energies. However, in the last few years the determinations of 
Aahad from (Q) have relied more and more on theoretical input. First perturbative 
QCD was used to better describe R(s) (for ^/s > 3 GeV) in energy regions above 
the resonances to the next flavour threshold Then, encouraged by the success of 
perturbative QCD to describe r decay, it was used down to 1.8 GeV, across a region 
with sparse data on R(s), giving || 

Aa had (M|) = (277.5 ± 1.7) x 10- 4 , (5) 

where ±10 comes from the uncertainty in perturbative QCD. 

Recent methods to calculate a(M|) have used the analytic behaviour of II (s) in 
the complex s plane. II is the hadronic contribution to the photon vacuum polarisation 
amplitude (or two-point correlator), where 

Aa had (s) = -47raReII(s), R(s) = 127rlmn(s). (6) 

Groote et al. H introduce special polynomial weight functions, which are arranged 
to suppress data contributions in specific regions and to essentially replace them by 
perturbative QCD (together with non-perturbative gluon and light quark condensate 
contributions). For example, for a specific interval (s a ,S(,), the idea is to deform 
the contour of integration so that the QCD contribution is evaluated around circular 
contours of radii s a and s&, away from singularities on the s a , Sb cut[|. The accuracy 
of this approach was improved in || to give an error on Aahad(-^l) °f ±1-6 x 10 -4 , 
close to the error of the straightforward estimate of (g) . 

We will see below that QCD does not reproduce the structure of the data for R(s) 
in the 1.8 < y/s < 3 GeV region. This region itself gives at least an error of ±10 -4 , 
which is additional to the errors of ||, ||. If the QCD description fails to reproduce 
the data, in an average sense in this interval, then nothing is gained by deforming the 
contour of integration^ using the trick of || g|. It is just not possible to circumvent 
these measurements of i?(s), and their uncertainties, in this way. 

In an interesting development Jegerlehner et al. [To| calculate the value of 
A«had at the large negative scale s — —M\ and then analytically continue around the 
large circle to positive s = M%. The QCD error comes mainly from the uncertainty 
in as(M^), and the neglect of a>g and higher order corrections. The extra error due 
to the analytic continuation around the large circle of radius ~ Mz is very small. At 
first sight it appears that we have dispensed with the data for R(s) altogether, and 
hence determine Aahad(Aff) with an error < 10~ 4 arising just from QCD. 

Unfortunately there are sizeable errors in II(— s) due to the uncertainty in what 
to take for the light quark masses; even the uncertainty in the charm threshold, Am 2 c , 
gives a sizeable error. Moreover the constant term in the oig contribution is not known 
yet. However these problems disappear for the derivative, the Adler function, 

D(-s) = -U^s^l (7) 
as 

f The QCD contribution was also evaluated around a circular contour in the s plane in Rcf. Jt|. 
X Due to Cauchy's theorem the result is independent of the choice of contour. 
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or for the discontinuity 

f)7T 

R(s) = — discll(s) = 12ImIT(s). (8) 

i 

Indeed this was recognized in || where D was evaluated for negative s using 
perturbative QCD (and parton condensate contributions). These authors checked 
that the theoretical values for D(— s) in the space-like region agree well with the 
"direct" evaluation in which the data for R(s) are used to determine A«had(— s). The 
agreement was found to be good even for y/s as low as 2 GeV. Unfortunately this is 
for the derivative, and not the value, of Aahad- To determine a(M|) Jergerlehner 
therefore evaluates 

Aa had (-M 2 z ) = [Aa had (-M|)-A« had (- So )] QCD + Aa had (- So ) data (9) 

where the QCD contribution is accurately known in terms of D and so the error in 
crfMz) dominantly reflects, once again, the error in the data for R(s). The evaluation 
@ can be done for any y/s~o > 2 GeV. Jergerlehner jlo| chooses y/so — 2.5 GeV and 
finds 

Aa had (Mf) = Aa had (-Afi) + (0.45 ±0.02) x 10~ 4 

= (277.82 ±2.54) x 10^ 4 , (10) 

where the error is entirely attributed to that on Aai iad (~so)- 

Since the error in a(M§ ) is driven by the data for R(s), it is relevant to discuss 
the new and forthcoming measurements. First, the CND-2 and SND experiments 
at VEPP-2M at Novosibirsk have made precision measurements of a wide range of 
exclusive processes in the interval 0.4 < y/s < 1.4 GeV @ @- 

For example, the 

measurement of e + e~ — > tt + it~ JTl| ] has considerably improved precision such that, 
together with chiral perturbation theory, we find that the error in the 2n contribution 
to Aahad for y/s < 0.96 GeV is reduced to about ±0.5 x 10 -4 , see Table 1. 

The second source of new information is the Beijing e + e" collider, where BES II 
are making direct measurements of R(s) for 2 < y/s < 5 GeV. So far values of R(s) at 
six values of y/s have been published [^3|. The two lowest values (y/s — 2.6, 3.2 GeV) 
are shown in Fig. 1, together with preliminary BES II data |H. These data, taken 
together with the compatible earlier 772 data [|l5|, indicate the existence of a "dip- 
bump" type structure in the measurements of R(s) about the 0(a^) perturbative 
QCD prediction shown by the continuous curve, see Fig. 1. The bump appears to be 
at y/s ^2.6 GeV and the dip at y/s ~ 2 GeV. This apparent structure leads to an 
uncertainty in Aah ad - Interestingly, the contributions to the right-hand-side of (^) 
from the interval 1.8 < y/s < 3 GeV obtained using for R(s) (i) the interpolation 
through the data and (ii) the perturbative QCD prediction are 

Aa ha 8 cT 3 GcV (Mz) = 17.77 ± 1.80 and 17.35 ±0.13, (11) 

respectively, in units of 10~ 4 . That is the results are very similar, which looks 
reasonable from the hadron-parton duality and QCD sum rule viewpoint. We take 
the data value and so, in fact, the uncertainty from this domain may be less than our 
pessimistic estimate. 

In the energy intervals described by perturbative QCD (3 < y/s < 3.74 GeV 
and y/s > 5 GeV) we use both the two-loop expression with the quark mass 
explicity included and the massless three-loop expression 111] calculated in the MS 
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Figure 1. A plot of the ratio R(s) versus y/s in the most sensitive ^fs interval. Up to 
\/5 = 2.125 GeV we show by continuous lines the upper and lower bounds of the sum 
of exclusive channels. Above y/s = 1.4 GeV we show the exclusive measurements 
of R(s), together with the interpolation (dashed curve) used in (4). We show six 
preliminary BES-II points, which were not used in the analysis, together with two 
published points which were used. The central perturbative QCD prediction is also 
shown. In the region 1.46 < \fs < 2.125 GeV we used, in turn, the inclusive and 



exclusive data to evaluate Act 



(■ r >) 



renormalization scheme|f|. In addition to the above we include the perturbative QCD 
error coming from varying m c ,mb, Mz within the uncertainties quoted in the |Tc| , 
as(M|) = 0.119 ± 0.002 and varying the scale as(cs) in the range 0.25 < c < 4§. 

f The uncertainty due to using a different scheme may be estimated to be of the order of the 
C(«g) correction, which is about 3 ^ e^r^^g /n) 4 . We may take r$ = —128 Jl7| which leads to an 
uncertainty much smaller than that given in Tabic 1. 

f We thank Thomas Teubner for valuable discussions concerning the perturbative QCD contribution. 
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By far the largest uncertainty comes from the region 1.4 < y/s < 3 GeV. The 
recent Novosibirsk measurements have considerably improved the exclusive data below 
1.4 GeV, but above 1.4 GeV we see from Fig. 1 that the sum of the exclusive channels 
exceeds the inclusive measurements of R. We therefore evaluated the contribution in 
the region 1.46 < y/s < 2.125 GeV using first the inclusive data| for R and then using 
the sum of the exclusive channels. (We checked that our results for the contributions 
from the exclusive channels were in agreement with the detailed table of results given 
in Ref. Jljj], if we were to omit the new Novosibirsk data and to include the r data.) 

Table 1. Contributions to Aoi la( j(s) X 10 4 of (^) coming from the different \fp 
intervals for three different values of s. The alternative values in the round brackets 
take the summation of exclusive channels as the contribution from the region 1.46 - 
2.125 GeV, rather than that from the inclusive data for R(s'). 



\fp interval 
(GeV) 


s = Ml 


s = -6 GeV 2 


a = -Afg 


2m* - 1.46° 
1.46 - 2.125 

2.125 - 3 
3 - 3.74 d 
3.74 - 5 

5 - oo d 
uj, (j), ^'s, T's 


38 - 18± {o°6 5 £ } 
11.80 ± 0.94 c 
(14.59 ± 1.76 b ) 
13.05 ± 1.22 c 
7.41 ± 0.04 
15.11 ± 0.50 
170.06 ± 0.66 
18.57 ± 0.57 


™ - { 25 } 

7.80 ± 0.60 c 
(9.59 ± 1.14 b ) 
6.33 ± 0.58 c 
2.58 ± 0.01 
3.64 ± 0.12 
6.09 ± 0.03 
11.26 ± 0.27 


38 - i8± {o :S } 

11.79 ± 0.94 c 
(14.58 ± 1.76 b ) 

13.03 ± 1.22 c 
7.39 ± 0.04 

15.04 ± 0.50 
169.83 ± 0.44 
18.52 ± 0.57 


Aa^d x 10 4 


274.18 ± 2.52 
(276.97 ±2.90) 


71.94 ± 1.40 
(73.73 ± 1.82) 


273.78 ± 2.47 
(276.57 ±2.85) 



a The upper (lower) error corresponds to the 2-7T (remaining) exclusive channels. 

b,c,d g rr0 rs with identical superscripts are added linearly. The remaining errors are added in quadrature. 



The results for A«had(s) of (|4|) are shown in Table 1 not only for s = M§, but 
also for two space-like values s = — so = — 6 GeV 2 and s = — M| in order to study 
Jegerlehner's approach |t0|. We see that indeed the error on the space- like evaluation 
at s = —6 GeV 2 is reduced in comparison to that for s = ±M%; as expected from 
the form of (^) the error mainly arises from uncertainties in the data for R(s') with 
s' < | so |- However before we can take advantage of this reduction we must consider 
the error in the perturbative continuation from s = — sq to — M§. (It is reasonable to 
assume that the error in going round the circle to s = M§ is much smaller.) To 
estimate the error in the difference Aa na d( — Mz) ~ Aa n ad(~ sq) we use the pure 
perturbative expression for R(s') in the whole interval of s' . The resulting uncertainty 
is the order of ±1 x 10~ 4 and so in practice there is only a marginal gain from the space- 
like evaluation. Until this theoretical error is more accurately quantified, we take a 
conservative approach and use the error of the s = Ml determination of Aa n ad(Ajf§ ). 

Let us return to the dip-bump structure in R(s) in the region 2 < \fs < 3 GeV, 
see Fig. 1. This was not expected, at least from the quark states; recall the success 

f The inclusive measurements of the 772 collaboration (with radiative corrections included) at higher 
energies appear consistent with the new Beijing measurements, which would tend to favour the 
inclusive input leading to a smaller value of Aa^ad ■ 



5 



of the perturbative QCD description of r hadronic decay. Charm is too heavy, and 
the light quarks too light to produce such a broad structure. Perturbative QCD was 
expected to reproduce the R(s) data in this y/s interval. However at C(a|) we open 
up the 3-gluon contribution, see Fig. 2, and we may expect some structure due to this 
channel. Indeed lattice computations predict a J p = glueball with mass typically 
in the range 3.3-4 GeV Q. So non-perturbative gluonic structure at C(a|), and 
higher order, may be anticipated in the region up to 4 GeV. In this connection we 
note that the non-perturbative (condensate) contribution to the Adler function D(—s) 
diverges for values of < 2 GeV, even in the space-like region M . Recall in the 





Figure 2. O(ctg) three gluon contribution to -R. The contribution of the light 
quarks cancel as e q = for q = u,d, s, so we are left with the c quark contribution. 



phcnomcnological description of the decay J/ip — > ggj that it was crucial to introduce 
an effective gluon mass of about 800 MeV to reproduce the 7 energy spectrum f2l|| . 
A similar result holds for T and Z decays [^2], [23]. We thus anticipate an effective 
threshold for the 3 gluon contribution to R(s) at y/s ~ 2.4 GeV. This effective gluon 
mass is also consistent with the predicted 1.6 GeV mass of the lightest (0 ++ ) glueball. 
More or less the same position of the threshold in R(s) is obtained assuming that the 
hadronization of the 3 gluon contribution to R (cut A in Fig. 2) proceeds through the 
6-quark channel. In the threshold region we must use the constituent quark mass, so 
that the 3 gluon channel opens at y/s > 6M cons tit — 2.1 GeV. 

Fig. 2 has not only the possibility of a gluon cut, but instead may be cut across the 
qq state, cut B. With respect to this latter cut, the imaginary part in the gluon channel 
plays the role of an absorptive correction. Therefore the diagram may give a negative 
qq contribution and a positive 3<? contribution leading to a non-trivial structure^. 
In addition there may be mixing, and mass shifts, of the J p — 1 glueball and qq 
resonant states. Of course the amplitude of this gluonic-driven structure cannot be too 
large as it originates from and higher order graphs. However the non-perturbative 
3 gluon interactions are strong enough to produce a resonant type structure in the 
relevant y/s region with up to about 10%| oscillations, and seen in Fig. 1. Note 
that the non-perturbative contribution to the Addler D function becomes large at 
y/\s\ ~ 2 GeV even in the space- like region |§. If the bump is due to 3-gluon 
production then we expect a larger yield of rf or 77 mesons (which are strongly coupled 
to gluons) in this energy region j24|. 

We conclude that it is important to use the experimental measurements of R(s) to 
calculate Aa na d and its error. One can replace the data for R(s) in a specific region by 

f Recall the cusp structure for a resonance occurring near a threshold in multichannel processes. 
I Estimated from the difference between O(a^) and 0{a 2 s ) QCD expectations. 
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the QCD prediction only after we understand the origin of all structure in the region, 
and are confident that our QCD approximation is sufficient to account for the effects. 
For example, much accuracy is gained by replacing the data measurements for R(s) 
above the T resonance region by the 0{a s ) QCD prediction. However we have seen 
that a similar replacement around the charm threshold is not so clear. The theoretical 
uncertainty due to the charm mass and the (more than two-loop) non-perturbative 
interaction may be larger than the accuracy of the data. From the presently available 
data we find (see Table 1) 

Aa had( M i) = 274 - 18 ± 2 - 52 or 276 - 97 ± 2 - 90 ( 12 ) 

according to whether the inclusive or exclusive data are used in the interval 1.46 < 
yfs < 2.125 GeV. When used in (Q) these values correspond to a value of the QED 
coupling at the Z pole of 

aiMlY 1 = 128.973 ± 0.035 or 128.934 ± 0.040 (13) 

respectively. It may be argued that the first value is favoured, since the 772 inclusive 
data appear to agree with the new Beijing measurements at higher energies, see Fig. 1. 
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